Magnetic phases of bosons with synthetic spin-orbit coupling in optical lattices 
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We investigate magnetic properties in the superfluid and Mott-insulating states of two- component 
bosons with spin-orbit (SO) coupling in 2D square optical lattices. The spin-independent hopping 
integral t and SO coupled one A are fitted from band structure calculations in the continuum, which 
exhibit oscillations as increasing the SO coupling strength. The magnetic superexchange model is 
derived in the Mott-insulating state with one-particle per-site, characterized by the Dzyaloshinsky- 
Moriya (DM) interaction. In the limit of |A| <C \t\, we find a spin spiral Mott state whose pitch value 
is the same as that in the incommensurate superfluid state, while in the opposite limit \t\ <C |A|, the 
ground state exhibits with a 2 x 2 in-plane spin pattern. 
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Quantum many-body states with spontaneous incom- 
mensurate modulated structures have attracted con- 
siderable interests in the past decades, and occur in 
many settings of condensed matter and ultracold atom 
physics. Celebrated examples include the incommen- 
surate magnetism with long-range and short range 
magnetic orders^, the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) pairing states 3,4 . Recently, the Bose-Einstein 
condensations (BEC) with spin-orbit (SO) coupling in- 
troduce a new member to this family. The SO coupled 
BECs are genuinely new phenomena due to the fact that 
the kinetic energy is not just a Laplacian but also linearly 
depends on momentum, which gives rise to the complex- 
valued and even quateronic valued condensate wavefunc- 
tions beyond the "no-node" theorem^. 

An important property of SO coupled condensates of 
bosons is that they can spontaneously break time-reversal 
symmetry, which is impossible in conventional BECs of 
both superfluid 4 He and many experiments of ultra-cold 
alkali bosons 7 . For example, it is predicted that such con- 
densates can spontaneously develop half-quantum vortex 
coexisting with 2D skyrmion-type spin textures in the 
harmonic trap£. Experimentally, spin textures of the SO 
coupled bosons have been observed in exciton conden- 
sations, which is a solid state boson systems with rel- 
ativist ic SO coupling 9 . Theoretically, extensive studies 
have been performed for SO coupled bosons which exhibit 
various spin orderings and textures arising from competi- 
tions among SO coupling, interaction, and confining trap 
energy^ 2 - - — . 

In the optical lattice, the SO coupled bosons are even 
more interesting. Early investigations have showed that 
the characteristic incommensurate wavevectors are in- 
commensurate with the lattice^. In this article, we study 
the SO coupled Bose-Hubbard model, focusing on the 
magnetic properties. The tight binding model is con- 
structed and the spin-independent hopping integral t and 
SO coupled hopping integral A are calculated as functions 
of the SO coupling strength in the continuum. Mag- 
netic superexchange models are derived characterized by 
the Dzyaloshinsky-Moriya (DM) interactio n 26 ! 27 . In the 



Mott-insulating phase, the single particle condensation 
is suppressed but the spin orderings are not. The spin 
orderings are solved in two different limits, |A| <C and 
|t| <C |A|, respectively. In the former case, the DM term 
destabilizes the ferromagnetic state to spin spirals, while 
in the latter case, a 2 x 2 in-plane spin ordering is formed. 

We begin with the non-interacting Hamiltonian of 
bosons with the Rashba SO coupling in a square lattice 
optical potential as 
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'-(ak x a y + Pk y a x ) + V(x,y), (1) 



where k so is the magnitude of wavevectors of laser 
beams generating SO coupling, a and f3 character- 
ize the anisotropy of SO coupling. Below we consider 
two situations. First, SO coupling is only along the in- 
direction, i.e., a = 1, /3 = 0, which agrees with the recent 
experiments^. Second, the isotropic Rashba SO cou- 
pling with a = 1, P = 1. V(x, y) is the periodic potential 
produced by laser beams with wavelength Aq as 



V(x, y) = — Vq cos 2 kox + cos 2 koy 



where fco = 27t/Aq, and the recoil energy E r = 



h z k z 
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We 



define a dimensionless parameter 70 = k so /ko to charac- 
terize the strength of SO coupling. The lattice constant 
& — Ao/2, and the reciprocal lattice is Gi = (^,0), 
G2 = (0, ^). The band structure of Eq. [I] is calculated 
by using the plane- wave basis. 

In the absence of SO coupling, the two-component 
bosons with strong optical potentials can be described 
by the lattice Bose-Hubbard model as 



H HU b = - Yl ^ [ 6 l A* + h A + H 
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where a =t, I denote the pseudospin components; 
and b\ a are bosonic annihilation and creation operators 
for spin a at site z, respectively. Yluj) denotes the sum- 
mation over all the nearest neighbors, rii is the boson 
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density operator at site i: rii — So-^LAo-- Generally, 
the interaction can be spin-dependent, and we only con- 
sider the spin-independent interaction below. Previous 
study of magnetism of two-component bosons includes 
Ref. [25], and we will consider it in the presence of SO 
coupling. 

The case of a = 1, j3 = is directly related to current 
experiments in the absence of optical lattice 11 . The SO 
coupling induces an extra term in the tight-binding term 
as 

H so = -A ^[^^4 - blju+g.d + hx, (4) 

i 

where e x is the unit vector along the x-direction. In 
momentum space, Eq. [4] becomes H so = X^k ^k^k^k? 

where \I/k = [&k,t> ^k,J T , anc ^ is a 2 by 2 matrix reads 
as: 

H\c = £\J + 2Asin/u £C (J 2/ (5) 

where = —2(t x cos k x + t y cos k y ) — fi and t x (t y ) is the 
hopping integrals along x and ^/-directions, respectively. 
In the long- wave limit k —> 0, Eq. (j5j) reduces to the 
Hamiltonian in continuous space realized in experiments. 

The SO coupling is equivalent to a pure gauge at /3 = 0, 
which can be eliminated by a gauge transformation 

U = exp{ik so • ra z }, (6) 

which applies to the doublet (6^,6^). Physically, this 
gauge transformation corresponds to a local pseudo-spin 
rotation by 0(r) = k so • r about the pseudo-spin z axis. 
The energy spectra of Eq. [5] has two branches as E± = 
—2t y [cos(k x ±k so )+cos k y ]—fjL, and the following relations 
are satisfied 

tx = ty cos k SOl A = ty sin k so . (7) 

Bosons condense into the energy minima of ±Q = 
(±fci,0) with ki = arctan(A/t 2/ ). The corresponding sin- 
gle particle wavefunctions at these two minima are: 

At the Hartree-Fock level, bosons can take either of \I/q sc 
as a plane-wave spin-polarized state, or, a superposition 
of them as -^(\I/q + \I/_q) = [cosQ • r, sinQ • r] T with 
the same energy. The latter one can be stabilized by 
spin-dependent interaction of H sPi i nt = U' XX n *,T — n hi) 2 
with U' < 0. It exhibits a spin spiral states in the xz- 
plane with the pitch wavevector 2Q as plotted in Fig. [TJ 
We will see that in the Mott-insulating state, although 
strong interaction suppresses the superfluidity, the spin 
configuration remains the same spiral order. 

We consider the Mott insulating state at (m) = 1, 
and construct the superexchange Hamiltonian for the 
pseudospin-^ bosons as 

H eff = J2[Hi, i+ e v +H i - l+ a (9) 



FIG. 1: Spin spiral configurations of the Bose-Hubbard model 
with unidirectional SO coupling. It is valid for both the in- 
commensurate superfluid state, and the Mott insulating state. 

For the vertical bond without SO coupling, Hi^+ ey is just 
the SU(2) ferromagnetic Heisenberg superexchang o 28 i 29 
as #i,i+e y = -Ji,ySi • S i+ e y where J ljV = U 2 y /U > 
0. For the horizontal bond, the SO coupling leads 
to the Dzyaloshinsky-Moriya (DM) type superexchange 
terms 26,27 as 

Hi,i+e x ~ —J\,x&\ ' Si+e K — ^12^^+^ • (Si X Si+g^ ) 

+ J 2 [Si • S i+ g x - 2(Si • d^i+^XSi+e, • di 5 i + ^)], 

(10) 

where J 2 = 4A 2 //7, J 12 = 4t x \/U. di, i+ e x is a 3D DM 
vector defined on the bond [i, i + e x ], and, di 5 i+e x = e y . 

The DM term of Eq.([TD|) prefers a spin spiral order- 
ing along the horizontal direction, as illustrated in FigJTJ 
The effect of the gauge transformation Eq. [6] on spin 
operators is to rotate Si around 7/-axis at the angle of 
2m6 where m is the horizontal coordinate of site i and 
= arctan(A/^)— , such that 

S'i = (l-cos2m(9)[d-Si]d + cos2m(9Si 

- sin2ra<9 Si x d, (11) 

where d = d{^+g x = e y . Through this transformation, 
the DM interaction is gauged away, and Eq. ([T0]) turns 
into a ferromagnetic coupling: 

#i,i+e x = -JoS[ • S- + ^, (12) 

where Jo = J\, x — 4(t 2 +A 2 )//7. The exchange model be- 
comes an isotropic ferromagnetic Heisenberg model, and 
thus spin polarization can point along any direction. In 
order to obtain the actual spin spiral configuration, we 
need to do the inverse operation of Eq. [TTJ Say, if we 
choose the classic spin at the point of origin along z di- 
rection S[ ,o] = &z-> according to the rotation defined in 
Eq. ([TT]h all the spins in the classic ground state are re- 
stricted within the x-z plane, and the classic spin at the 
point [m, n] is S^ m ^ = cos(2m0)e z + sm(2m0)e x . As 
shown in FigJT] (b), the classic spins form a chiral pat- 
tern with a characteristic length, which is the same as in 
the superfluid case as plotted in Fig. [TJ The only differ- 
ence is that the superfluid phase coherence is lost in the 
Mott-insulating state. 

Now we discuss the isotropic Rashba SO coupling with 
a = f) = 1. From the symmetry analysis, we easily have 
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FIG. 2: A) The dependence of the spin-independent hopping 
integral t and the spin-dependent one A v.s. the SO coupling 
strength 7 = k so /ko. The optical potential depth is Vb = SE r . 
B) Sketch of Wannier wavefunctions for f(r) (solid black line), 
f(r — a) (solid red line) and g(r) (dashed blue line) in Eq. [16] 
where a is the lattice constant 

t x = t y for spin independent hoppings, while the spin- 
dependent SO hoppings become 

h'so = -*J2{ b U b i+^- b L b i+^] +h - c > 

- iX J2l b ti b i+^ + bl t h +gy J + h.c. (13) 

i 

In momentum space, the tight-binding band Hamiltonian 
turns to: H' = ^ k ^k#k^k, where 

ff k = e^I + 2X[sm k x a y + sin k y a x ], (14) 

where = — 2t(cosk x x + cos k y y). The energy spectra 
of Eq.Q read 



E± = e k ± 2Ay sin 2 ^ + sin 2 fc y . (15) 

The band minima are four-fold degenerate at the points 
Q sc = (±fc, ±fc), where fc = arctan 

Next we calculate the band parameters t and A versus 
SO coupling parameter 7, by fitting the band spectra us- 
ing the plane-wave basis in the continuum. The results 
are plotted in Fig. [2] A. Both t and A oscillate and de- 
cay as increasing 7, which can be understood from the 
behavior of the onsite Wannier functions. Each optical 
site can be viewed as a local harmonic potential and the 
lowest single particle state wavefunction was calculated 
inRef. @ 

^=i(r) = [/(r)^(r)e^] T , (16) 

and its time-reversal partner is i/jj z= _i(f) = 

(—g(r)e l< P,f(r)). f(r) and g(r) are real radial wavefunc- 
tions, which exhibit characteristic oscillations with the 
pitch value k so and a relative phase shift approximately 
I as plotted in Fig. [2] B. t and A are related to the 
off-centered integrals of f(r) and g(r) of two sites, 
which overlap in the middle. As a result, t and A also 
oscillate as increasing 7, which also exhibit a phase shift 
approximately at ir/2 as shown in Fig. [2] A. 



We would like to clarify one important and subtle 
point. Actually the on-site Wannier functions are no- 
longer spin eigenstates, but total angular momentum 
eigenstate j z = |, and thus are still a pair of Kramer 
doublets. For the operators (&^,6^) T defined on site 
z, they do not refer to spin eigenbasis but to the j z - 
eigenbasis. In fact, in the case that k so > ^o, the onsite 
spin moments are nearly zero. The j^-movements mainly 
come from orbital angular momentum. As pointed out 
in Ref. [8[, the Wannier functions of j z eigenstates ex- 
hibit skyrmion-type spin texture distributions and half- 
quantum vortex on each site. This phenomena also re- 
mind us of the Friedel oscillation in solid state physics. 
In the case of k so ^> &o, each site exhibits Landau level- 
type quantization: states with different values of j z are 
nearly degenerated^, and a single band picture ceases to 
work here. 



(a) (b) 




FIG. 3: (A) The pattern of DM vectors of the superexchange 
magnetic model in the Mott-insulating state. (B) Illustration 
of the frustration in the spin configuration with DM interac- 
tions, the rotations around x-axis and ^/-axis do not commu- 
tate with each other. 

Deep inside the Mott-insulating phase, we obtain the 
effective magnetic Hamiltonian: 

H[ i+ £ is the same as Eq. (fT0|) . and 

H l,i+e y = -JlSi ' Si+e y - Jl2di,i + e y ■ (Si X S i+ gJ 

+ J 2 [Si • S i+ e, - 2(Si • di, i+ eJ(S i+ e, ■ 4,1+^)], (18) 

where (1^+^ = & x . The pattern of the DM vectors is 
shown in Fig|3] (a), which is a strongly reminiscent of 
that in cuprate superconductor YBa2Cu306 31,32 . The 
classical ground state of Eq. ([TT|) is nontrivial because the 
DM interaction can not be gauged away: DM vectors in 
horizontal bonds favor spiraling around the i/-axis, while 
that in vertical bonds favor spiraling around the x-axis. 
Since rotations around x and y-axis do not commune, 
no spin configurations can simultaneously satisfy both 
requirements, which leads to spin frustrations shown in 
Fig. 0(B). 

The quantum situation of Eq. [TT] is even more in- 
volved, which can only be solved approximately. Below 
we focus on two situations: |A| <C |t| and |A| ^> At 
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A = 0, the ground state of Eq. JT7J) is known to be 
ferromagnetism. At \X/t\ <C 1, we use spin- wave approx- 
imation to analyze the instability of a ferromagnetic state 
induced by the DM interaction. Notice that in this case, 
it is impossible to find a global rotation as in Eq. ([TT]) 
to gauge away the DM vectors and transform Eq. (fTTj) 
to an SO (3) invariant Hamiltonian, thus the quantized 
axis in the spin wave analysis can not be chosen arbitrar- 
ily. To gain some insight, we choose a classic ferromag- 
netic state as a variational ground state parameterized 
by Si = S (cos 'y sin 77, sin sin 77, cos 77) . The correspond- 
ing variational energy Eq = —S 2 (J\ — J2 + 2 sin 2 r]J 2 ) 
is minimized when 77 = 7r/2, which implies that the xy- 
plane is the easy plane. 
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FIG. 4: (a) Spin spiral ordering in the limit of |A| <C \t\. (b) 
The 2x2 pattern in limit of \t\ < |A|. 

To calculate the spin waves spectra, it is convenient to 
rotate the coordinate so that the new z-axis points along 
the direction 1 = [110] in the original coordinate (we 
choose 1 as the quantized axis). The Holstein-Primakoff 
transformation is employed to transform Eq. (fT7|) into 
the bosonic Hamiltonian: 

H b = ^-J (cos20 - ism20/V2)a\ai+ e ^ + h.c, (19) 

Z,/i 

where = arctan(A/t) as defined above, \i = x, y. We 
only keep quadric terms and ignore the terms propor- 
tional to sin 2 since X/t <C 1. In momentum space, it 
becomes 

H' ex = —2 Jo [cos 20 cos k x H — — sin 20 sin 
k V 2 

+ cos 20 cos k y H — ^= sin20sinfc 2/ ] cj^Ck, (20) 
v2 



The minimum of the dispersion of Eq. [20] occurs at 
points Qm = (±k f ,±k f ), where k' satisfies that tan A: 7 = 
-j= tan 20. Compare it with the energy minima in the 
noninteracting band Hamiltonian Q s / = (±fc, ±fc), we 
have k' — 2k at the limit of 7 — » 0. The nonzero mini- 
mum of the magnon spectrum is a signature of the spin 
spiral order, as shown in Fig. HI (a). 

In the opposite limit of \X/t\ ^> 1, if we consider the 
Eq. (fTTj) as a classic Hamiltonian, the J2 term has a 
continuous degenerate manifold. The DM J12 energy is 
zero for all of them, thus does not lift the degeneracy. The 
ferromagnetic coupling J\ term selects the configuration 
of spin lying the x?/-plane with the 2x2 pattern shown 
in Fig. Hl(b). This is a discrete symmetry breaking state, 
thus there is no Goldstone mode for this 2x2 states. The 
magnon excitation should be gapped even at quantum 
level, and such a state is stable at large values of X/t. 
For the intermediate values of X/t, we expect quantum 
phase transitions from the spin spiral state in Fig. HI (a) 
to the 2x2 state in (b). Due to the very different spin 
configurations, rich phase structures with exotic patterns 
of spin spirals and textures are expected. 

The magnetic orders proposed above can be detected 
by the spin-dependent Bragg scattering. The transition 
rate for the elastic Bragg scattering is directly related to 
the spin structure factor for the atoms, thus will exhibit 
peak in the characteristic momentum Q of the spin spiral 
states. Similar methods have been proposed to detect the 
antiferromagnetic stated as well as the FFLO stated in 
cold atoms. 

In conclusion, we have investigated the magnetic or- 
dering of two-component Bose-Hubbard model with syn- 
thetic SO coupling. The band parameters of hopping in- 
tegrals exhibit characteristic oscillations as increasing SO 
coupling strength, and the onsite magnetic moments are 
nearly orbital moments at large SO coupling strength. In 
the Mott-insulating state with one particle per site, an ef- 
fective magnetic superexchange model with the DM type 
interaction is derive. The spin spiral state and the 2x2 
states are found in the limits of |A| <C \t\ and |A| ^> \t\, 
respectively. 
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Note added Upon the posting of this paper, we become 
aware two papers on the similar topi o 35 i 36 . After the 
paper was posted, two papers also appear 37 
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